Energy-resolved spatial inhomogeneity of disordered Mott systems 
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We investigate the effects of weak to moderate disorder on the T = Mott metal-insulator 
transition in two dimensions. Our model calculations demonstrate that the electronic states close to 
the Fermi energy become more spatially homogeneous in the critical region. Remarkably, the higher 
energy states show the opposite behavior: they display enhanced spatial inhomogeneity precisely in 
the close vicinity to the Mott transition. We suggest that such energy-resolved disorder screening 
is a generic property of disordered Mott systems. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.23.-k, 71.30.+h 



I. INTRODUCTION 

The observation of a metal-insulator transition in high- 
mobility two-dimensional electron systems at zero mag- 
netic field has sparked renewed interest in this type of 
transition^. In these systems, electron-electron interac- 
tions represent the largest energy scale in the problem^. 
Further evidence for the crucial role of electronic cor- 
relations has come from complementary experiments^ 
reporting a substantial mass enhancement close to the 
metal-insulator transition. Taken together these ex- 
perimental results stress the importance of the ill- 
understood effects of disorder in strongly correlated elec- 
tronic systems^. 

In one of the first studies of its kind, Tanaskovic et 
al. have investigated the interplay of strong correla- 
tions and disorder using a dynamical mean field the- 
ory (DMFT)^ approach^. The DMFT approach to dis- 
ordered systems treats correlations on a local level by 
solving the embedded-atom strongly correlated Ander- 
son impurity problem in the self-consistently determined 
fixed bath of the other electrons. One is thus forced to 
consider an ensemble of single-impurity actions, one for 
each lattice site. Remarkably, Tanaskovic et al. found 
that very strong site disorder screening emerges precisely 
in the vicinity of the Mott metal-insulator transition^. 
This effect can be traced back to the pinning of the single- 
impurity Kondo resonances to the Fermi level, which acts 
to suppress the effective randomness. 

Motivated by this striking result, we extended their 
work to finite dimensions. To treat both strong cor- 
relations and disorder in a non-perturbative fashion at 
T = 0, we use a generaHzation of the DMFT, the statis- 
tical DMFT (statDMFT)2., implemented using a slave 
boson impurity solver—. This approach is equivalent 
to the description of the effects of disorder through a 
Gutzwiller-type wave functio n^°i^^ , which in the clean 
case realizes the Brinkman-Rice scenario of the Mott 
transitions^. The statDMFT method retains the local 
treatment of electronic correlations. However, in con- 
trast to the infinite-dimensional DMFT approach, here 
each strongly correlated site sees a different hath of elec- 



trons, refiecting the strong spatial fiuctuations of their 
immediate environment. We find that, concomitant to a 
strong pinning effect, additional effective disorder is also 
generated through an increasingly broader distribution 
of quasiparticle weights. Their interplay leads to a non- 
trivial landscape in energy space: the proximity to the 
Mott transition acts to suppress density of states fiuctu- 
ations close to the Fermi level, while at the same time 
enhancing them at higher energies. 



II. STRONGLY CORRELATED THEORY 



We focus on the disordered Hubbard model 



H ^ e.,c\^c^a - t ^ (cLcja + h.( 



(1) 



where c|^(cio-) is the creation (annihilation) operator of 
an electron with spin projection a on site z, —t is the 
nearest-neighbor hopping amplitude, U is the on-site 
Hubbard repulsion, nia = c\^Cia is the number operator, 
and the site energies are uniformly distributed in the 
interval [— W/2, VF/2], where W is the disorder strength. 
We work at half fiUing (chemical potential /i = U/2) on 
an LXL square lattice with periodic boundary conditions. 
All energies will be expressed in units of the clean Fermi 
energy (the half-bandwidth D) Ep = it. 

We treat the Hamiltonian of Eq.[T]in its paramagnetic 
phase within the statDMFT—. This theory is exact in the 
non-interacting Hmit and reduces to the standard DMFT 
in the absence of disorder. UnHke the DMFT, however, 
it incorporates Anderson localization effects. We start by 
writing an effective action (in imaginary time) for a given 
site i, with the simplification that we neglect all non- 
quadratic terms in the local fermionic operators except 
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for the local /7-term 
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The site i is connected with the rest of the lattice through 
the bath (or "cavity") function (r) , which in stat- 
DMFT (but in contrast to DMFT) varies from site to 
site and thus exhibits strong spatial fluctuations. 

The effective action in Eq. [2] is precisely the action 
of an Anderson impurity modeU^ embedded in a sea 
of conduction electrons described by (t). Therefore, 
this approach maps the original Hubbard Hamiltonian 
in Eq. [T] onto an ensemble of single-impurity Anderson 
Hamiltoniansi^. The local i-site Green's function, calcu- 
lated under the dynamics dictated by the effective action 
in Eq. [21 can be written as {iuj is a Matsubara frequency) 



El {iuj) ~ A, {iu>) ' 



(3) 



which also serves as a definition of Si (iuj) , the z-site self- 
energy. It is important to point out that within stat- 
DMFT the electronic self-energy (iuj) is still local, al- 
beit site-dependent. 

The bath function A^ (iuj) can be viewed as the Weiss 
field of this mean field theory, here elevated to a full 
function of frequency or time. It is determined through a 
self-consistency condition that demands that the Green's 
function G-°^ (iuj) obtained from the effective action in 
Eq. [2] be equal to the diagonal (local) part of the full 
lattice Green's function 



Gii [iuj) 
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e - Ho - S (iw) 



(4) 



where S [iu) and e are site-diagonal matrices 
\T,{iuj)\j = Ei(iw)(5y, [e]j^- = StSij and Ho is the 
clean {W = 0) and non-interacting {U — 0) lattice 
Hamiltonian. In general, this step involves the inversion 
of the frequency-dependent matrix within brackets in 
Eq. Igl). 

It is worthwhile to point out that the statDMFT ap- 
proach requires a massive numerical effort since it is nec- 
essary to solve a single impurity problem for every lattice 
site as well to perform the inversion implied by the self- 
consistency condition. On the other hand, it provides 
access to entire distribution functions and accounts for 
spatial correlations between local quantities. 

To find Si {iuj), we need to solve the auxiliary single 
impurity problems for a given set of A^ {iuj). For this 
task, we have used the four-boson mean-field theory of 
Kotliar and Ruckenstein2. at T = 0, which is equivalent to 
the well-known Gutzwiller variational approximation. In 



practice, we need to solve a pair of non-linear equations 
for the site-dependent Kotliar-Ruckenstein slave boson 
amplitudes and dS^ 



^i^G^n^-) - Z.{l-el + d^, (5) 
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Eqs. ((MSI) involve the local i-site Green's function 
G-°'^ (icj). In the KotHar-Ruckenstein theory, the z-site 
self-energy Si (zcj) is given by 



Si (icj) = (l - Z- ^) + - ei + ^, 
which, when plugged into Eq. ^ yields 

Z, 



«w - ZiVi - ZiAi (iuj) ' 



(9) 



(10) 



It is clear that Zi has the physical interpretation of a 
quasi-particle weight (wave function renormaHzation) . It 
also renormaHzes the hybridization function, thus setting 
the local Kondo temperatur o^^i^^i^^ . Furthermore, from 
Eq. ([HI we see that Vi can be viewed as a renormalized 
on-site disorder potential 



+ S, (0) - fi. 



(11) 



We numerically solved the disordered Hubbard model 
with StatDMFT, using the KotHar-Ruckenstein theory as 
the impurity solver, for several lattice sizes up to L = 50. 
For every (C/, W) pair we typically generated around forty 
realizations of disorder. We carefully verified that for 
such large lattices, all our results are robust and essen- 
tially independent of the system size (see, e.g., the inset 
of Fig. [31). 



III. THE DISORDERED MOTT TRANSITION 

According to the scaling theory of localization, any 
amount of disorder drives a system with dimension equal 
to or smaller than two to an insulating phasei^. How- 
ever, this result was obtained in the absence of electron- 
electron interactions. In the last few years, strong nu- 
merical evidence has been obtained indicating that in- 
teractions can act to enhance the conducting proper- 
ties of two-dimensional electronic system o^^'^"'^^ . More- 
over, the stability of a 2d metal with respect to weak- 
localization corrections has been investigated by Pun- 
noose and Finkelsteir^ in very careful recent work. 
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Figure 1: Strength of the renormalized site energy disorder, as 
given by the standard deviation a (u) of P {v), normalized by 
its non-interacting value (see text for definition of v). Close 
to the Mott transition, the disorder screening remains strong 
even for moderate values of disorder. Results are shown for 
L = 20. In the inset we show the typical (Ztyp) and average 
(Zav) values of the local quasiparticle weight Zi as functions 
of the interaction U. The Mott transition is identified by the 
vanishing of Ztyp. We note that Zav is finite at Uc indicating 
that a fraction of the sites remains nearly empty or doubly 
occupied. Results are shown for L = 20 and W = 2.25. 



These authors have demonstrated that any 2d metal re- 
mains stable with respect to sufficiently weak disorder 
due to additional (anti-localizing) interaction corrections. 
Note, however, that both the weak-localization and the 
corresponding interaction corrections are manifested only 
through a very weak, logarithmic dependence on the sys- 
tem size. Such subtle finite size effects are not visible 
for the very weak (renormaHzed) disorder we deal with 
in this work. Nevertheless, based on the very convincing 
considerations of Punnoose and Finkelstein, the stability 
of a 2d metal is beyond immediate doubt and these issues 
are not of relevance for the questions we focus on in this 
paper. 

The low energy behavior on the metaUic side of the 
transition is characterized by the distribution of the lo- 
cal quasiparticle weight Zi. Following previous studies of 
the disordered Mott transition^i^^, we choose the typical 
value of Zi , here defined through the geometrical average 
Ztyp — exp {(InZi)}, as the order parameter of the tran- 
sition. This quantity vanishes linearly at a critical value 
of interaction Uc = Uc (W) at which the Mott transition 
takes place, marking the transmutation of most itinerant 
electrons into local magnetic moments (see the inset of 
Fig. [T|). Because random site energies tend to push the 
local occupation away from half fiUing, Uc{W) is an in- 
creasing function of the disorder strength W^. We also 
stress that the average value of the quasiparticle weight 
Zav is small yet finite at the Mott transition (inset of 
Fig. [U , indicating that some sites remain either empty 
(e^ = 1) ou doubly occupied (d^ = 1) and do not give rise 
to localized magnetic moments. 



IV. PARTICLE-HOLE SYMMETRY AND 
STRONG DISORDER SCREENING 

Using the fact that in the current statDMFT approach 
the lattice problem is mapped onto an ensemble of auxil- 
iary Anderson impurity problems, we can characterize 
the approach to the Mott transition by a steady de- 
crease of the local Kondo temperature (T^ oc Zi). 
Moreover, the renormalized disorder potential Vi can be 
thought of as the position of the local Kondo resonance 
energy. 

In the DMFT limit {d oo), each site has many neigh- 
bors and thus is embedded in the same (self-averaged) 
environment described by Aav (i^)- In this regime, as 
long as Aav (w) is particle-hole symmetric, we find per- 
fect disorder screening close to the Mott transition^. This 
happens because, as we increase U towards UdW), we 
approach the Kondo limit, Zi 0, causing Vi (the Kondo 
resonance) to be "pinned" to the Fermi energy^ 

However, within the statDMFT approach, (w) fiuc- 
tuates strongly from site to site and is not locally particle- 
hole symmetric (see Fig. [2]). For this reason, we expect 
that Vi will also have a contribution which is propor- 
tional to Re[Ai(0)]. Therefore, we have no guarantee 
that a similar mechanism of disorder screening will per- 
sist close to the critical point in 2d. Surprisingly, how- 
ever, for small and moderate disorder {W < Uc), we do 
get a very strong disorder screening close to the Mott 
metal- insulator transition (see Figs. [1] and [3]). 

V. MOTT DROPLETS 

The site-to-site fiuctuations in the hybridization func- 
tion encode spatial correlations between the physical 
quantities, implying, for example, that the value of 
the local quasi-particle weight at a given site depends 
on the quasi-particle weight values at the neighboring 
sites. Therefore, due to rare disorder configurations. 
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Figure 2: Real (dashed lines) and imaginary (full lines) parts 
of the hybridization function Ai (uj). Although the average 
value of Ai (u) is particle-hole symmetric, for an arbitrary 
site Ai (uj) is locally particle-hole asymmetric. Results are 
shown for L = 20 and W = (7 = 2.25. 
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Figure 3: Probability distribution function P (u) of the renor- 
malized site energy v. As we move towards the Mott transi- 
tion, disorder screening takes place and P (v) becomes in- 
creasingly narrower. Results are shown for L = 50 and 
W = 2.25. The inset illustrates how for such large lattices 
our results for P (v) are essentially independent of the system 
size: results are shown for L = 20, full line, and L = 50, 
dashed line at = 2.25 and U/Uc = 0.93. 



it is possible to find regions containing sites in which 
Zi <C Ztyp (red regions in Fig. |4]). Since the approach 
to the Mott insulator corresponds to Z ^ 0, such re- 
gions with Zi <^ Ztyp should be recognized as "almost 
localized" Mott droplets (consisting of localized magnetic 
moment o^^i^^i^^ ) inside the strongly correlated metallic 
host. Within our Brinkman-Rice picture, each local re- 
gion provides^ii^ a contribution Xi ~ 7j ~ ^o the 
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Figure 4: Spatial distribution of the of the local susceptibility 
Xi ~ normalized by its typical value Xtyp- From bot- 

tom to the top we have U/Uo = 0.73, 0.87, 0.96. The color 
scale is logarithmic in order to stress that, as we approach 
the critical point, we have the formation of regions in which 
Xi ^ Xtyp (localized magnetic moments) and Xi ^ Xtyp (An- 
derson insulator droplets). Results are shown for W = 2.25 
and L = 50. 
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Figure 5: Distribution of the quasiparticle weight divided by 
its typical value for W = 2.25 and (W) ~ 3.7. As we 
approach the critical point, P (Z/Ztyp) becomes increasingly 
broader, even though the bare disorder strength W is kept 
fixed. Results are shown for L = 20. 



spin susceptibility or the Sommerfeld coefficient, respec- 
tively. The local regions with the smallest Zi thus dom- 
inate the thermodynamic response and will ultimately 
give rise to an Electronic Griffiths Phase in the vicinity 
of the disordered Mott metal-insulator transitio n^^i^^i^° . 

The Mott droplets have a direct influence on the 
probability distribution function P (Z/Ztyp), since they 
give rise to a \ow-Z tail, contributing to the fact that 
P{Z/Ztyp) actually broadens as we approach the Mott 
metal- insulator transition for a fixed W, as shown in Figs. 
HandU 

Surely, there is also a contribution to the broadening 
of P{Z/Ztyp) coming from the high-Z tail originating 
from those sites with Zi ^ Ztyp. Such sites retain a 
finite Zi at the transition and are nearly empty or dou- 
bly occupied, giving rise to Anderson insulating regions 
(blue regions in Fig. The coexistence of locaHzed 
magnetic moments, and nearly empty or doubly occu- 
pied sites close to the critical point is characteristic of a 
two-fiuid behavior—. Nevertheless, this interesting phe- 
nomenon is not relevant for the present analysis, since 
the major contribution to the energy resolved inhomo- 
geneities, discussed in the next Section, as well as to the 
thermodynamic response, comes from the Mott droplets. 



VI. ENERGY-RESOLVED INHOMOGENEITIES 

From the above discussion, it is clear that the behaviors 
of P {v) and P{Z/Ztyp) near the Mott metal-insulator 
transition are quite distinct. While the former exhibits 
strong disorder screening (Fig. [31, the latter suggests 
that the disorder is actually increasing (Fig. This 
dichotomy gives rise to an energy- dependent effective dis- 
order^ which manifests itself, for example, in the spatial 
structure of the local density of states. The local density 
of states is defined as Pi (w) = (I/tt) Im [Ga (w — iO+)], 
where the lattice Green's function was given in Eq. lU 
and has the following expression within our Brinkman- 
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Rice picture 



Gii {lo) 



Z 



(12) 



where Z and v are site-diagonal matrices such that 
[Z]^- = ZiSij, [v],y = Vt6ij. Therefore, the frequency- 
dependent effective disorder potential "seen" by the quasi- 
particles at energy uj can be defined as 



(uj) ^ V,- — 



(13) 




LO* = uj/Ztyp, which we will keep constant as we approach 
the transition. In Fig. [6l we show topographic maps of 
the local density of states for one specific realization of 
disorder. Because of strong disorder screening, w 
close to the critical point. Thus, if the system is exam- 
ined at the Fermi energy {uj = 0), it becomes more and 
more homogeneous as the transition is approached. At 
higher energies {uj* ^ 0), however, the fiuctuations in 
Ztyp/Zi come into play. Since they are very pronounced 
close to the Mott transition, we instead find a strong 
enhancement of the spatial inhomogeneity. 

This result is surprisingly reminiscent of recent spec- 
troscopic images on doped cuprates^. Our theory, which 
focuses on local (Kondo-like) effects of strong correla- 
tions (while neglecting inter-site magnetic correlations) 
and does not include any physics associated with super- 
conducting pairing, strongly suggests that such energy- 
resolved inhomogeneity is a robust and general feature of 
disordered Mott systems. Indeed, recent results corrobo- 
rate this picture of an energy-resolved strong-correlation 
driven disorder screening^^. 
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Figure 6: Spatial distribution of the local density of states nor- 
malized by its clean and non-interacting value for one given 
realization of disorder and two distinct values of interaction: 
(a) U/Uc = 0.87 and (h) U/Uc = 0.96, where Uc ~ 3.7. We 
have w = in the bottom figures and uj/Ztyp = 0.1 in the 
top ones. At the Fermi energy, the local density of states 
distribution becomes homogeneous as we approach the Mott 
transition. Conversely, if we move even slightly away from the 
Fermi energy, the distribution becomes in fact more inhomo- 
geneous close to Uc- Results are shown for L — 50, W — 2.25. 

In our Brinkman-Rice scenario of the Mott transi- 
tion, the quasi-particle bandwidth is reduced by Ztyp as 
U Uc- To monitor the behavior within the quasi- 
particle band, we therefore introduce a rescaled frequency 



VII. CONCLUSIONS 

We have discussed the results of a Brinkman-Rice ap- 
proach to the disordered Mott transition. A striking fea- 
ture that is apparent in this scenario is the different be- 
haviors of the effective site disorder and the quasipar- 
ticle weights as the transition is approached. Whereas 
randomness in the former is suppressed, the latter be- 
comes extremely singular and broad on the way to the 
Mott insulator. The end result of this dichotomy is a 
non-trivial energy-space landscape, which is signaled by 
an energy-dependent inhomogeneity. 
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